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We investigate chiral symmetry breaking and strong CP violation effects in the phase diagram 
of strongly interacting matter. We demonstrate the effect of strong CP violating terms on the 
phase structure at finite temperature and densities in a 3-flavor Nambu-Jona-Lasinio (NJL) model 
including the Kobayashi-Maskawa-t'Hooft (KMT) determinant term. This is investigated using an 
explicit structure for the ground state in terms of quark-antiquark condensates for both in the scalar 
and the pseudoscalar channels. CP restoring transition with temperature at zero baryon density is 
found to be a second order transition at 6 — -k while the same at finite chemical potential and small 
temperature turns out to be a first order transition. Within the model, the tri-critical point turns 
out to be (T c , fj, c ) ~ (273, 94) MeV at 6 = -n for such a transition. 

PACS numbers: 12.38.Mh, 12.39.-x, 11.30.Rd, 11.30.Er 



I. INTRODUCTION 

Strong interaction is known to respect space and time reflection symmetry to a very high degree. However this is 
not a direct consequence of laws of quantum chromodynamics (QCD), which, in principle permit a parity violating 
term or the so called 0-term given as 
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; g 2 F"F a ^. (1) 



In the above, F* is the gluon field strength and F^ v being its dual. This term while being consistent with 
Lorentz invariance and gauge invariance, it violates charge conjugation and parity unless 9 = mod n. However, 
CP symmetry conserving nature of QCD has been established by precise experiments that sets limit on the intrinsic 
electric dipole moment of neutron. The current experimental limit on this leads to limit on the coefficient of the CP 
violating term of the QCD Lagrangian as < 0.7 x 10 -11 [l|. This smallness of the CP violation term or its complete 
absence is not understood completely though a possible explanation is given in terms of spontaneous breaking of a 
new symmetry : the Peccei-Quinn symmetry [2] which could give rise to axions. For zero temperature and and zero 
density, spontaneous parity violation does not arise for 9 = by the well known Vafa-Witten theorem Q. On the 
other hand for 9 = n there could be spontaneous CP violation by the so called Dashen phenomena [J- Because of 
the nonperturbative nature of this 9 term in QCD, this problem has been studied extensively in low energy effective 
theories like chiral perturbation theory linear sigma model @ as well as Nambu-Jona-Lasinio (NJL) model and 
its different extensions 

Even if CP is not violated for QCD vacuum, it is possible that it can be violated for QCD matter at finite 
temperature or density. It has been proposed that hot matter produced in heavy ion collision experiments can give 
rise to domains of meta stable states that violate CP [111 ], Experimental signatures for the existence of local CP 
violation has been based on charge separation of hadronic matter due to the strong magnetic field produced in heavy 
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ion collision experiments by a mechanism called chiral magnetic effect (CME) [12j . This mechanism may explain the 
charge separation in the recent STAR results [13( . 

In the present work we focus our attention on how chiral transition is affected when there is CP violating term in the 
Lagrangian. For this purpose, we adopt the 3-flavor NJL model as an effective theory for chiral symmetry breaking 
in strong interaction [3, [la]- The CP violating parameter 6 is included in the Kobayashi-Maskawa-t'Hooft (KMT) 
determinant term. In this context we note that the two flavor scenario for spontaneous CP violation for 8 — ir has 
been studied in this model |8fl . This has been further extended to study the restoration of CP at finite temperature 
The effect of the theta vacuum on the deconfinement and chiral transition has also been analyzed within a two 
flavor NJL model with Polyakov loop • 

We organize the present work as follows. In the next section we shall discuss the 3-flavor NJL model with a CP 
violating term. We consider a variational ground state with quark-antiquark pairs that is related to chiral symmetry 
breaking. The ansatz functions are to be determined through minimization of the thermodynamic potential. The 
ansatz is general enough to include both scalar as well as pseudoscaler condensates. As we shall see the pseudoscalar 
condensates develop for non zero values of 8 in the KMT determinant term. In section III we discuss the resulting 
phase diagram at finite temperature as well as finite density for different values of the CP violating parameter in the 
Lagrangian. In Section IV we summarize our results and give a possible outlook. 



II. NJL MODEL WITH CP VIOLATION AND AN ANSATZ FOR THE GROUND STATE 

To describe the chiral phase structure of strong interactions including the CP violating effects, we use the 3-flavor 
NJL model along with the flavor mixing determinant term. The Lagrangian is given by 

8 

£ = -m)ip + Gj2 [(^V) 2 + (^7 5 A A V) 2 ] - K [e i8 dei{$(l + 7 5 )V'} + e~ ie det{$(l - 7 5 )V>}] , (2) 

,4=0 

where tp l ' a denotes a quark field with color 'a' (a = r, g 1 b), and flavor '£' (i — u, d, s), indices. The matrix of current 
quark masses is given by m — diag/(m u , mj, m s ) in the flavor space. We shall assume in the present investigation, 
isospin symmetry with m u = rrid- In Eq.@, A" 4 , A = 1, • • • 8 denote the Cell-Mann matrices acting in the flavor space 

and A = ^/| If, where 1/ is the unit matrix in the flavor space. The four point interaction term ~ G is symmetric 

under SU(3)v x SU(3)a x U(1)v X U(1)a- The determinant term ~ K, which generates a six point interaction for 
the case of three flavors, breaks U(1)a symmetry for vanishing 8 values. The effect of topological term of Eq.([T]) is 
simulated by the determinant term of Eq.@ in the quark sector. 

We shall next consider an ansatz for the ground state with quark-antiquark condensates which includes both the 
scalar as well as CP violating pseudoscalar channel. To make notations clear, we first write down the field operator 
expansion for the quark fields as given in [l6l . Kt\ , 

lKx,t = 0) = J Vi(kK k - x dk = J [C/ (k)g°(k) + V (-k)g°(-k)] e^dk, (3) 

where C/°(k) and V°(— k) are the four component spinors which can be explicitly written as, 

™- (,.!£&>) ^-(^Sf)- <4) 

The superscript '0' indicates that the operators q° and q° are the two component operators for the quark annihilation 
and antiquark creation corresponding to the the perturbative or the chiral vacuum |0). Here we have suppressed the 
color and flavor indices of the quark field operators. The function x°(k) in the spinors in Eq.Q are given as 
cotXi* = Wi/|k|, for free massive fermion fields, i being the flavor index. For massless fields x°(|k|) = tt/2. 

We next consider an ansatz of the ground state at zero temperature as 



\n) = u q \o), 



(5) 
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where, U q — U q iU q ii is an unitary operator. U q i and U q n are unitary operators described in terms of quark- 
antiquark creation and annihilation operators. Explicitly they are given as 



U qI = exp Q rfk g °(k)+(o- • k) rs /(fc)^(-k) - h.c 



(6) 



and 



U qII = expfj dkq r (k)^rg{k)q- r (-k) - h.c^j (7) 

where f(k) and g(k) are the ansatz functions which are to be determined later from the extremization of the 
thermodynamic potential. 

Finally, to include the effect of temperature and baryon density, we use the techniques of thermofield dynamics 
(TFD) which is quite convenient while dealing with operators and states [3, [l9j]. Here, the statistical average of an 
operator is given as an expectation value over a 'thermal vacuum'. The methodology of TFD involves the doubling 
of the Hilbert space [ill]. Explicitly, the 'thermal vacuum' is constructed from the ground state at zero temperature 
and density through a thermal Bogoliubov transformation given as 

|fi(/0,M)) =U F \Sl) =e B(0^-B(P^)\n) (8) 

with, 

B\p,n) = [ [Mk, J 8,/iV(k)V(-I0 t +e+(k,p, f i)q'(k)q'(-k)]dk. (9) 



In Eq.(|9]) the ansatz functions 9±(k, j3, fx) will be related to the quark and the antiquark thermal distributions 
respectively and the underlined operators are the operators in the extended Hilbert space associated with thermal 
doubling in TFD method. 

In the following section we shall compute the thermodynamic potential which will involve calculating thermal 
average of different operators which is given by the expectation values of the corresponding operator with respect to 
the state given in Eq. (JEJ)- This can be evaluated directly by realizing that the state |0(/3, //)) is obtained from the 
state |0) by successive Bogoliubov transformation. Thus e.g. we have 

(0(/3 M )|Vl ia (x)<(y)|)fi(^)) = ^5 ab J ^ e -^(x-y) Al(kj/3j/i)/3a (10) 

where, 



A*(k,^) = i [(cos 2 ^ + sm 2 9 l _) + (sin 2 6> 4 1 - sin 2 6^X7° cos ft cos 2g l + a • ksin^ cos 2g l - i 7 ° 7 5 sin2 5 ,; (11) 



where, we have introduced a new function <fi l (k) = Xo + 2/ 4 (k) in terms of the condensate function pik) of Eq.©. 
From Eq. ([T(])) . it is easy to calculate the scalar and pseudoscalar condensates. In terms of the ansatz functions <fi l {k) 
and g l (k) the scalar and pseudoscalar condensates for the i-th flavor can be respectively written as 

2 AT f 

= J dkoos^' cos 2g\l -n l _ - n\) = -P s (12) 

. 2N 



2N 

(^75^), = -ijKZ^ I ^ksin2 3 4 (l - n l _ - n\) = -il l p (13) 



where n 1 ^ — sin 2 0V . Thus a non vanishing P s will imply chiral symmetry breaking phase while a non vanishing 
Ip or equivalently g l (k) will indicate CP violating phase. The condensate functions pik), g l (k) as well as the 
thermal functions ^(k, /3,/x) shall be determined by extremization of the thermodynamic potential with respect to 
the respective functions. We shall carry out these extremization in the following section. 
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III. EVALUATION OF THERMODYNAMIC POTENTIAL AND GAP EQUATIONS 

As mentioned we shall be considering the chiral phase structure in presence of the CP violating terms within the 
frame work of Nambu-Jona-Lasinio model of Eq.([2|). The energy density is given by the expectation value of the 
Hamiltonian corresponding to the Lagrangian given in Eq.@ with respect to the thermal ansatz of Eq.©. The 
energy density can be written as 



e = T + V = T + V s + V D 
where T is the expectation value of the kinetic term in Eq.(|2]) and using Eq. (|T0|) . is given as 



(14) 



2N f 



(15) 



Vs the contribution from the four point interaction term in Eq.([2]l to the energy density and using Eq. (|10[) , this is 
nven as 



(16) 



A=0 



Finally, Vd denotes the contribution from the determinant interaction term in Eq.® and using Eq. (|10p . is given as 



V D = K{e t0 detail + 75)^} + e^ 9 det{i>{l - 75)^}) 



2K 



cos < — 



2 \ e ijk \I s IpIp 



I P T3 
t ijk\ 1 s 1 s 1 p 



Now, the thermodynamic potential is given as 



£1 = e — pp 



S 



where p is the quark chemical potential corresponding to the quark number density p given as 



(17) 



(18) 



P 



Finally, S is the entropy density given as 



= E (^ = t|% E /dk(l-sin 2 ^+sin 2 ^). 



l — U.d.S 



(19) 



2N f 

S = dk (cos 2 9t In cos 2 6>1 + sin 2 6 l _ In sin 2 6 l _ + cos 2 6\ In cos 2 9\ + sin 2 6\ In sin 2 0+) . (20) 

Thus the thermodynamic potential given in Eq. (fT8|) is known in terms of the ansatz functions of Eq.®. Extremizing 
the thermodynamic potential with respect to </> l (k) and g l (k) respectively leads to 



tan < 



K 

Mi 



and tan 2g l 



ML 



ML 1 + lkl 2 



(21) 



where M* and M % v are respectively the contributions to the constituent quark mass (for i-th flavor) from the scalar 
and pseudoscalar condensates and they are given by, 
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Ml = m* + AGI\ + K \e l3k \{cos9{IlI k s - iffi) - sin6{P s I k p + P p I k s )} (22) 
M; = 4G7;-if|e^|{cos^/ p fc +7|/ s fc )- S m^/ p fe -7^)}. (23) 

Finally extremizing the thermodynamic potential with respect to the thermal distribution functions leads to 



sin 2 9\ = — , (24) 

± exp(w l Tfi') + 1 v ' 



where, w l (k) = Jk 2 + (Mf + Mf). Thus we can see that the constituent quark masses get contribution from 
both the scalar and pseudoscalar condensates. 

Substituting the extremized solution for the condensate functions tan</>* and tan2g I in Eq. (fT3|) and Eq. ([T3|) . we 
have the self consistent equations for the scalar and pseudoscalar condensates, 



- IN f A/P 

11 S = (# J dk ^ ~ H - ^ (25) 

in r M i 

P p ee 1(4,^4,), = y dk (1 - nl - n\) -f . (26) 

Thus with the scalar and the pseudoscalar condensates given as above, Eq.s (|23H23f are actually coupled self- 
consistent equations for M] and M p . 

Substituting the extremized solutions for the condensate functions and using the gap equations Eq. (|23H23"|) in 
Eq. dlSp . the thermodynamic potential becomes 



9. 



dktu* - w) + 2 Gs y [if - if] + e K 1 ; 



(27T) 

3 

+4K cos 6 'JJ/j - 27fsin( 



i=i 



3 i 



i=l 



2A^ C 
>(2tt) 



dk [lnll + e-^^'-^^ + lnjl + e 



-/3(^'+p l )-i 



(27) 



In the above we have subtracted the perturbative vacuum energy density contribution. Let us note that, the effective 
potential has been calculated using an explicit ansatz for the condensate and not evaluating the effective potential 
at a mean field level after performing a chiral transformation for the quarks so as to remove it from the determinant 
term as it has been computed in Ref. [1, [13] ■ 

Eq. (|27j) for the thermodynamic potential and the gap equations for the scalar and pseudoscalar masses i.e. Eq. 
[2^| shall be the focus of our numerical analysis which we do in the following section. 



IV. RESULTS AND DISCUSSIONS 



For numerical calculations, we have taken the values of the parameters of the NJL model as follows. The coupling 
constant G s has the dimension of [Mass]~ 2 while the six fermion coupling K has a dimension [Mass] -5 . To regularize 
the divergent integrals we use a sharp cut-off, A in 3-momentum space. Thus we have five parameters in total, namely 
the current quark masses for the non strange and strange quarks, m q and m s , the two couplings G s , K and the 
three-momentum cutoff A. We have chosen here A = 0.6023 GeV, G S A 2 = 1.835, A'A 5 = 12.36, m q = 5.5 MeV and 
m s = 0.1407 GeV as has been used in Ref. [Til] . After choosing m q — 5.5 MeV, the remaining four parameters are fixed 
by fitting to the pion decay constant and the masses of pion, kaon and r/' . With this set of parameters the mass of rj 
is underestimated by about six percent and the constituent masses of the light quarks turn out to be M u,d = 0.368 
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Fig. 1-a Fig. 1-b 

FIG. 1: ^-dependence of the condensates (Fig. 1-a) and the effective potential at zero temperature and zero baryon density 
(Fig. 1-b). 

GeV for u-d quarks, while the same for strange quark turns out as M s = 0.549 GeV, at zero temperature and zero 
density. 

For a given temperature and the chemical potential, we first solve the coupled self consistent gap equations Eq.s (|231 
|2"3")) with the parameters of the model as above. Since we have assumed isospin symmetry and have m u = m d , these 
are actually four coupled equations: two for the scalar condensates related to the two masses M" = Mf , M| and 
two equations for the pseudoscalar condensate related to the corresponding mass parameters = M d , M*. The 
solutions to these equations are then substituted in Eq. (|27[) and checked regarding minimum of the thermodynamic 
potential. If there are more solutions to the gap equation, the one with the minimum thermodynamic potential is 
chosen. 

Let us first discuss the ground state structure at zero temperature and zero density. In Fig. (1-a) we show the theta 
dependence of contributions to the mass of up quark from the scalar as well as the pseudoscalar condensates. As 
is clearly seen as 9 increases the condensates in the two channels behave in a complimentary manner. While the 
magnitude of scalar condensates decreases with 9 (till 9 = n) , the magnitude of the pseudoscalar condensate increases 

so that the total constituent quark mass M = ^Jm s 2 + M p 2 remains almost the same. Spontaneous CP violation is 
clearly seen for 9 — it with two degenerate solutions for differing by a sign. In Fig. (1-b) we show the effective 
potential as calculated above as a function of 9. The effective potential is normalized with respect to the same at 
9 = 0. The minimum of the potential is at 9 = which is consistent with the Vafa-Witten theorem and has a cusp at 
9 = 7r which has also been observed in 2-flavor NJL model [8J. 

Next we consider the effect of nonzero density and temperature. In Fig.©, we show the variation of masses of the 
quarks with chemical potential at zero temperature. In Fig.(2-a) we show the variation of masses for the case 9 = 0. 
In this case the pseudoscalar condensates vanish and the contribution to the masses of the quarks are from the scalar 
condensates only. The (approximate) first order chiral transition takes place at \x ~ 361 MeV for u and d quarks 
with their masses decreasing discontinuously to about M u ' d ~ 52 MeV from their vacuum value of M u,d = 368 MeV. 
Because of the flavor mixing KMT term, this decrease is reflected also in the decrease of the strange quark mass to 
M s = 464 MeV from its vacuum value of M s = 549 MeV. This result is similar to the results obtained in the context 
of color superconductivity in NJL model with a determinant term (20| and in the context of chiral symmetry breaking 
in a similar model [21). Similarly in Fig.(2-b) we show the variation of masses of up and strange quark as well as 
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FIG. 2: Quark masses as a function of quark chemical potential at zero temperature for 8 — (Fig. 2-a), 9 — tt/2 (Fig. 2-b) and 
— it (Fig.2-c). The pseudoscalar contribution for strange quarks is zero in this range of chemical potential. 



the variation of the contributions from scalar and pseudoscalar condensates to the constituent mass of up quark for 
8 = tt/2. Here the critical chemical potential for chiral transition is /i c ~ 375 MeV where the mass contributions from 
scalar and pseudoscalar condensate become 24 MeV and 12 MeV respectively from their vacuum values of 266 MeV 
and 248 MeV. The total mass for the u and d quarks become 27 MeV from its vacuum value of 364 MeV. On the 
other hand the contribution to the strange mass from the pseudoscalar condensate is negligible (~ 12 MeV) compared 
to the contribution from the scalar condensate (~ 548 MeV). Because of flavor mixing again, strange quark mass also 
decrease to 463 MeV at fi c — 375 MeV. For 8 — tt the scalar condensate almost vanishes but for the nonzero current 
quark masses while the contribution to the constituent quark mass arises from the pseudoscalar condensate as shown 
in Fig. (2-c). As the quark chemical potential is increased there is a first order transition at [i c ~ 368 MeV. At [i c the 
pseudoscalar condensate vanishes and the contribution to quark mass arises solely from the scalar condensate which 
is non vanishing because of the nonzero current quark masses. 

Next we discuss the condensate variations with temperature. For 8 = 0, at zero baryon density the chiral crossover 
transition takes place for temperature about 200 MeV as may be seen in Fig.(3-a). As 8 increases, the pseudoscalar 
condensates starts becoming nonzero and increase with 8. For 8 = tt/2, masses arising from both type of condensates 
are shown in Fig.(3-b). Here, both the scalar and pseudoscalar masses show a crossover transition as temperature 
is increased. In Fig.(3-c), we show the behavior of the masses for 8 = tt. The transition for the pseudoscalar mass 
becomes a second order transition at 8 = tt instead of a crossover which was the nature of transition for lower 
8. This feature is elaborated in Fig. (J^J) where 8 dependence of the nature of transition of scalar and pseudoscalar 
condensates with temperature is shown for zero baryon density. Fig.(4-a) shows the 8 dependence of the transitions 
for the scalar condensate. The transition is always a crossover for scalar condensate. Fig.(4-b) shows the transitions 
of the pseudoscalar condensates for different 8. We can see that the transition is a second order transition for 
8 = tt whereas it is a crossover for other values of 8. Similar kind of results have been obtained in sigma model 
calculations @ but there the transition at 8 = tt is a first order transition instead of a second order transition. 
The CP restoring transition temperature for zero baryon density case turns out to be 192 MeV. However, the total 
constituent mass is nonzero as the scalar condensate is non vanishing again because of nonzero current quark masses. 
This high temperature restoration of CP is expected as the instanton effects responsible for CP violating phase become 
suppressed exponentially at high temperature [22| . 

The CP restoring transition in the present 3-flavor NJL model turns to be second order for zero chemical potential 
case similar to the results for the 2-flavor case @ unlike the case of linear sigma model coupled to quarks [f| . The reason 
behind such different behavior regarding the order of the transition is due to the non analytic vacuum term in the NJL 
model [9j. However, a first order CP transition is observed with finite chemical potential and small temperature. This 
is clearly shown in Fig. ([5]) where we have shown the dynamical mass, arising from the pseudoscalar condensate, 
for different temperatures as a function of the quark chemical potential for 8 = tt. While at zero temperature, the 
order parameter decreases discontinuously, as the temperature increases, it becomes less sharp and finally results in 
a second order transition at high temperature. 
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FIG. 3: Quark masses as a function of temperature at zero quark chemical potential for 9 = (Fig.a), 6 = 7r/2(Fig b) and 

6 = TV. 
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FIG. 4: The nature of transition of the scalar (Fig.a) and pseudoscalar (Fig.b) condensates with temperature at zero baryon 
density for different values of 9. 



In Fig.©, we show the phase diagram in the plane of quark chemical potential and temperature for the CP violating 
transition. Since the transition is first order at zero temperature and second order at zero chemical potential, there 
is a tricritical point for this transition in this plane. This turns out to be (/j, c ,T c )— (273,94) MeV. Including Polyakov 
loop for the two flavor NJL model, such a tricrtical point occurs at (209,165) MeV [lfj. First order transitions are 
associated with existence of meta-stable states. CP is restored in these meta stable states and these are the nontrivial 
solutions of the gap equation Eg. (1231) . however with lower pressure than the stable solutions. In the phase diagram 
of Fig.©, such solution exist in the region between solid line and the dotted line. 
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FIG. 5: Contribution from the pseudoscalar condensates to the u-quark mass as a function of quark chemical potential for 
different temperatures. Here we have taken 8 — ir. 
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FIG. 6: The phase diagram for CP transition in the T — fi plane. The region between the solid line and the dotted line have 
solutions to the pseudoscalar mass gap equation but with higher thermodynamic potential. Here we have taken 8 — ir. 



V. SUMMARY 



To summarize, we have tried to examine the effect of 9 vacuum on the phase diagram of strong interaction. This is 
investigated within the framework of a 3-fl.avor NJL model. The effect of CP violating 9 term of QCD is incorporated 
through the KMT determinant interaction term in the quark space. 

In the two flavor NJL model, it has been observed that spontaneous CP violation for 9 = ir occurs depending on 
the magnitude of current quark masses as well as the strength of the determinant coupling 8] . In the present case of 
3-flavor NJL model, we observe that spontaneous CP violation occurs for 9 = tt for the phenomenologically consistent 
parameters [l5| for the current quark masses as well as the strength of the determinant coupling. 



10 



To calculate the thermodynamical potential in presence of a 9 term, instead of performing a chiral rotation of the 
quark field operators [13] , we have used a variational approach using an explicit construct for the ground state. We 
have considered an ansatz state general enough to have condensates both in scalar as well as pseudoscalar channel. 
The ansatz functions are determined through minimization of the thermodynamic potential. 

Apart from the temperature that has been considered earlier for two flavor cases [f| Q we have also considered the 
effect of finite quark chemical potential and discussed the phase diagram in the T — fi plane. It turns out that for 
the range of temperature and the chemical potentials we have considered, the strange quark condensates do not get 
dynamically generated in the CP violating pseudoscalar channel even for non vanishing 9. None the less the strange 
quark antiquark condensates in the scalar channel do affect the pseudoscalar light quark condensates through the 
flavor mixing coupled gap equations. 

The CP restoring transition turns out to be a crossover for non vanishing values of 9 that become a second order 
transition for 9 = ir for zero chemical potential. However, such a transition is a first order transition at small enough 
temperature and as the quark chemical potential increases. This leads to a tricritical point in the the phase diagram 
for CP transition. 
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